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ABSTRACT 

A  computer  model  was  developed  to  simulate  a  C.R.P.  propeller 
system.   A  set  of  simultaneous  differential  equations  was  solved 
using.  CSMP  III-the  IBM  simulation  language.   The  computer  model  was 
verified  against  manufacturer's  data.   Power  piston  displacement, 
power  piston  velocity  and  chamber  pressure  were  determined  for  a 
nominal  supply  pressure  of  600  PSI  and  various  loadings.   Overall, 
the  system  model  and  data  should  prove  useful  for  investigations  of 
design  variations. 
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NOMENCLATURE 

A  cross-section  area 

B  load  damping  coefficient 

C  leakage  coefficient 

F  force 

K  load  spring  constant 

K  valve  pressure  coefficient 
c  c 

K  valve  flow  gain 

q 

Kf  flow  force  spring  rate 

M  mass 

P  load  pressure  on  power  piston 
m  *  * 

P0  load  pressure  on  input  piston 

P  supply  pressure  on  power  piston 

P  supply  pressure  on  input  servovalve 

S3. 

0  flow  rate 

t  time 

V  volume 

w  valve  area  gradient 

X  input  piston  displacement 

X  valve  displacement 

Vcl 

Y  power  piston  displacement 

S  effective  bulk  modulus 
e 

C  damping  factor 

to,  natural  frequency 


SUBSCRIPTS 

a  refers  to  input  piston 

£  refers  to  load 

m  refers  to  power  piston 

p  refers  to  piston 

s  refers  to  supply  pressure 

v  refers  to  valve 
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I.   INTRODUCTION 

During  recent  years,  naval  applications  of  gas  turbine  engine 
systems  have  become  more  prominent.   At  first,  the  gas  turbine  engine 
was  developed  as  an  aircraft  power  source.   Presently,  it  is  entering 
service  as  a  marine  power  source.   Gas  turbine  engine  systems  have 
several  advantages  compared  to  other  power  systems.   The  main  advan- 
tages are  its  relative  compactness  and  the  fact  that  the  rotating  out- 
put of  the  turbine  is  usually  what  is  wanted  to  turn  a  wheel,  propeller, 
generator,  or  pump.   Unfortunately,  compared  to  other  systems,  there 
is  a  problem  of  matching  the  characteristics  of  the  load  and  the  engine. 
In  this  respect,  controllable  reversible  pitch  (C.R.P.)  propellers  are 
useful. 

Controllable  pitch  propellers,  unlike  fixed  pitch  propellers,  have 
the  ability  to  reverse  thrust  through  a  reversal  of  pitch  rather  than 
a  reversal  of  rotation.   The  controllable  pitch  propeller  is  much  more 
complex  than  a  solid  propeller  and  contains  many  more  parts  [1] .   How- 
ever, the  operation  of  a  C.R.P.  system  provides  a  greater  degree  of 
flexibility  in  the  maneuvering  and  docking  of  ships. 

Since  there  is  an  increasing  interest  in  the  application  of  con- 
trollable reversible  pitch  propellers  to  naval  vessels,  the  objectives 
of  this  study  were  as  follows: 

1.  to  determine  an  overall  analytical  model  of  an  operational 
C.R.P.  system; 

2.  to  check  the  system  response  for  design  variations. 


11 


II.   BACKGROUND 

A.   BASIC  CONTROL  CONCEPTS 

In  a  control  system  in  which  the  output  appears  as  mechanical  energy, 
the  quantity  to  be  controlled  can  be  either  load,  torque,  position,  or 
velocity  and  can  be  controlled  under  either  open-  or  closed-loop  condi- 
tions.  In  an  open-loop  system  [2] ,  the  input  is  fed  to  the  controller 
which  then  operates  on  the  output.   The  input,  being  an  independent  vari- 
able, is  incapable  of  correcting  any  errors  between  the  desired  and 
actual  output. 

The  accuracy  of  an  open-loop  system  is,  therefore,  greatly  affected 
by  external  disturbances  and  variations  of  system  parameters.   In  a  closed- 
loop  system  [2] ,  the  output  is  fed  back  and  compared  with  the  desired 
value.   The  error  between  the  desired  and  actual  output  is  fed  into  the 
control  elements  which  in  turn  operate  on  the  output,  tending  to  reduce 
or  null  the  error.   This  is  the  concept  of  negative  feedback.   A  block 
diagram  of  a  closed-loop  system  is  shown  in  Figure  1. 


REFERENCE 
INPUT 


ERROR 
SIGNAL 


MANIPULATED 
VARIABLE 


CONTROLLER 


ACTUATOR 


CONTROLLED 
VARIABLE 


FEEDBACK  SIGNAL 


FEEDBACK 
ELEMENT 


Figure  1.   Block  Diagram  of  a  Closed-Loop  System. 
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In  a  position  control  system  in  which  the  manipulated  variables  are 
pressure  and  flow,  the  hydraulic  quantities  are  converted  into  a  dis- 
placement by  action  on  the  piston  rod  or  output  shaft.   The  two  proper- 
ties of  flow  which  can  be  manipulated  are  its  direction  and  magnitude. 
Directional  control  of  flow  is  used  for  simple  position  control  systems 
under  both  open-  and  closed- loop  conditions  [3] . 


Return 


Figure  2.   Valve  Controlled  Piston. 


Directional  control  valves  that  could  be  used  here  are  of  the  three- 
or  four-way  type.   A  schematic  arrangement  of  a  four-way  valve  control- 
ling a  symmetrical  piston  is  shown  in  Figure  2. 

When  the  control  valve  as  shown  is  in  the  neutral  position,  the 
piston  is  hydraulically  locked.   A  displacement  of  the  valve  spool  in 
either  direction  causes  the  piston  to  move  in  the  opposite  sense. 
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An  idealized  graph  of  a  position  control  system  which  is  operated 
under  open-loop  conditions,  by  a  directional  control  valve  is  shown  in 
Fiaure  3. 


Figure  3.   Position  Control  System  Response  (Open-Loop) . 

When  the  valve  is  displaced  from  its  neutral  position,  the  actuator 
moves  at  constant  velocity  in  the  required  direction  until  the  valve 
is  returned  to  neutral.   An  open-loop  system  of  this  type  is  only  suit- 
able for  simple  applications  which  do  not  require  accurate  piston  control, 

Closed-loop  position  control  systems  employ  controllers  having  non- 
continuous  or  "on-off"  characteristics.  When  the  actuating  signal 
reaches  a  threshold  value,  the  directional  control  valve  becomes  opera- 
tive and  causes  the  actuator  to  move  at  constant  velocity  in  the  chosen 
direction.  It  continues  to  move  until  the  error  signal  falls  below  the 
threshold  value.  This  causes  the  valve  to  return  to  its  closed  position 
and  stop  the  motion  of  the  piston. 
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The  most  widely  used  type  of  control  for  hydraulic  positioning  servo 
systems  is  the  control  flow  type,  where  both  flow  direction  and  magnitude 
are  manipulated.   The  two  methods  of  achieving  flow  control  are: 

1.  controlling  the  displacement  of  a  variable  delivery  pump; 

2.  throttling  fluid  through  the  variable  orifices  of  a  control 
valve. 

The  two  basic  differences  between  pump  operated  and  valve  operated 
servo  systems  can  be  described  as  follows: 

1.  In  pump  operated  servos,  working  pressure  is  variable, 
depending  on  the  external  loading,  while  valve  controlled 
servos  operate  at  constant  system  pressure. 

2.  In  valve  controlled  servos,  power  sources  are  outside 
the  loop. 

In  this  study,  a  valve  controlled  servo  system  was  studied. 

B.   ANALYTICAL  APPROACH 

This  analytical  investigation  was  directed  toward  predicting  various 
C.R.P.  system  pressures  and  motions.   An  analytical  model  was  developed 
using  first^order  differential  equations.   The  mathematical  model,  which 
was  programmed  on  a  digital  computer,  predicted  the  position  and  velocity 
of  the  power  piston  and  the  pressure  at  various  points  in  the  system  as 
a  function  of  time.   Further  analysis  was  performed  using  the  model  to 
obtain  preliminary  predictions  of  system  performance  under  a  variety  of 
external  loadings. 

To  understand  the  derivation  in  the  following  sections,  Figure  4 
indicates  the  geometry  of  a  basic  controllable  reversible  pitch  (C.R.P.) 
propeller  system.   Beginning  with  the  input  side,  a  force  balance  on  the 
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input  servovalve  yielded  an  equation  for  the  load  pressure  on  the  input 
piston  in  terms  of  the  valve  motion.  The  input  piston  was  moved  by  load 
pressure,  P0  .   The  valve  rod  control  assembly  was  moved  by  the  force 

363. 

exerted  on  the  input  piston.   This  caused  the  primary  supply  pressure,  P  , 

to  pass  through  the  hub  valve  into  the  power  piston  and  move  the  propeller 

crank-slot  pin.   This  power  piston  control  system,  as  such,  is  similar  in 

operation  to  a  vehicular  power-assisted  steering  mechanism. 

To  facilitate  the  digital  simulation  it  was  assumed  that  the  input 

valve  spool  displacement,  X   ,  and  the  arbitrary  load,  F0 ,  were  step 

va  J6 

inputs.   In  the  process  of  solution,  the  system  state  equations  were 
solved  using  CSMP  III, the  IBM  simulation  language. 
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III.   CONSTITUENT  RELATIONSHIPS 

It  was  a  difficult  problem  to  identify  the  dynamic  characteristics 
of  the  controllable  pitch  propeller  system  because  of  its  essentially 
nonlinear  behavior.   Identification  was  partially  achieved  by  assuming 
that  these  responses  were  linear  and  then  applying  parameter  estimation 
for  linear  systems. 

A.   INPUT  SERVOVALVE  AND  PISTON 

The  combination  of  a  servovalve  and  a  piston  is  a  common  hydraulic 
power  element  [3] .   In  analyzing  this  simple  system,  three  basic  rela- 
tionships surface.   The  first  concerns  the  valve,  the  second  concerns 
the  piston  dynamics  and  the  third  concerns  the  flow  of  hydraulic  fluid 
between  the  elements. 

Regarding  the  valve,  the  basic  equation  that  pertains  is 


Q0   =  K   X   -  K   Pg  (1) 

36a    qa  va    ca  Jca 


where  K   =  70  w  J& 

qa       a   sa 


K 


=  ^a 

"ca   P 

sa 


Regarding  the  piston,  the  basic  equation  governing  value  dynamics  is 


ZF  =  m  X  (2) 

p    v 
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The  sum  of  the  force  acting  on  the  piston  is 


ZF      =   F    0    -   F0       -   K        X   -   B         X  (3) 

p  pic  £a  pa  pa 


where  F   „    =  A     P„ 

pJ6         p     ica 

F„       =   K       (X   -   Y)    =    0.86  w    (P      -   P    ) (X   -    Y) 
36a  r  m     s  m 


Combining  equations  (2)  and  (3)  gives 


A   Pn   =  M   X  +  B    X  +  K    X  +  F„  (4) 

p  £a    v      pa      pa      ica 


The  flow  relationship  between  the  two  elements  is  developed  by 
applying  the  continuity  equation  to  each  of  the  piston  chambers.   In 
final  form,  this  continuity  relation  becomes 

2£a  "  CtPJla  P*a  +  Ap  *  +  (7T»  *ta  (5) 

e 

where  C   »   is  an  overall  leakage  coefficient. 

Equations  (1) ,  (4)  and  (5)  are  the  basic  equations  of  this  sub- 
system. 

B.   POWER  PISTON 

From  the  characteristic  equation  of  the  critical  center  spool  valve 
(hub  valve) ,  one  can  derive 

Q_  =  K    (X  -  Y)  -  K   P  (6) 

Tn    qm  cm  m 

where       K   =  140  w   /p 
qm        m   s 

cm   P 
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In  a  manner  similar  to  the  development  of  equation  (5) ,  the  continuity 
expressions  for  the  hub  valve  can  be  combined  to  yield 


V 

0   =  C.  0      P   +A   Y  +  ("7"?-)  p  (7) 

~m    tp£m  mm      4  $    m 

where  Co   is  an  overall  leakage  coefficient. 

The  final  equation  arises  by  applying  Newton's  second  law  to  the 
power  piston  [4] .   The  resulting  force  equation  is 


A  P   =  M„  Y  +  B    (X  -  Y)  +  K    (X  -  Y)  +  Fn        (8) 

m  m    I  pm  pm  I 


Equations  (6) ,  (7)  and  (8)  are  the  basic  equations  of  the  power 
piston  system. 

C.   OVERALL  SYSTEM 

In  analyzing  the  overall  system,  two  modes  of  operation  became 

apparent.   First,  if  the  system  input  was  a  forward  valve  displacement, 

X   ,  the  output  was  a  power  piston  position  and  velocity.   The  second 

mode  occurred  during  recovery  from  an  external  loading.   For  this  mode, 

the  input  was  an  external  load,  FQ .   At  this  point,  it  was  necessary 

to  simplify  the  governing  equations.   Usually,  the  loading  on  a  valve 

controlled  piston  system  is  simpler  than  the  most  general  in  two 

respects  [3] .   First,  significant  spring  loads  are  usually  absent  and 

hence  K   =  K   =0.   Second,  damping  loads  are  likewise  negligible  and 

then  B   =  B   =0. 
pa    pm 
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Thus,  equation  (8)  reduced  to 


*  =  ~TT     {Am     Pm  "  F?>  {9) 

M„    m  m    x, 


A  combination  of  equations  (6)  and  (7)  produced 

4  6 

P   =  — (-K     P-AY  +  K    X-K    Y)         (10) 

m    V       cem  m    m      qm      qm 
m 

Equation  (4)  was  simplified  to 


5  =  -5"  (-  K£  X  +  Kf  Y  +  Ap  PJ.a>  (11) 

V 


Finally  equations  (1)  and  (5)  yielded 


4  $ 

p0  =  —rr^   (_K   p„-ax  +  k   x)         (12) 

la.  V       cea  £a    p      qa  va 


Equations  (9-12)  were  the  governing  equations  for  the  C.R.P.  system. 
The  following  is  a  listing  of  the  constants  used  in  the  derivation  of 
these  governing  equations  and  in  the  accompanying  computer  solution. 
The  valves  were  derived  from  manufacturer's  data  [5]. 
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TABLE  I 


LIST  OF  CONSTANTS 


A 
P 

(AA) 

A 
m 

(AM) 

K       (CEM) 
cm 

K 

ca 

(CEA) 

L 

V 
a 

(VA) 

V 
n 

(VM) 

w 
a 

( ) 

w 

m 

(  W) 

P 
s 

(PS) 

P 

sa 

(PSA) 

6e 

(BE) 

\ 

(XM) 

M 
V 

(VM) 

K 

qa 

(QAK) 

t 

K       (QMK) 
qm   x 

cross-sectional  area  of  input  piston 

cross-sectional  area  of  power  piston 

flew  gain  of  input  servovalve 

flow  gain  of  power  piston 

total  volume  of  input  piston 

total  volume  of  power  piston 

valve  gradient  of  input  servovalve 

valve  gradient  of  power  piston 

supply  pressure  (power  piston) 

supply  pressure  (input  servovalve) 

bulk  modulus 

total  mass  of  power  piston  and 
propeller 

total  mass  of  input  piston  and 
valve  rod  assembly 

total  flow  pressure  coefficient  of 
input  piston  and  servovalve 

total  flow  pressure  coefficient  of 
power  piston 


57.67  in 


1102.0  in 


1.54  in3/sec/PSI 

3 
0.11  in  /sec/PSI 


749.71  in' 


14326.0  in 
1.571  in2/in 


11.938  in  /in 
400.0  PSI 
600.0  PSI 
100000.0  PSI 
69.8757  lb-sec2/in 

2 

4.9  lb-sec  /in 

2057.08  in3/sec/in 

3 
33426.54  in  /sec/in 
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IV.   SYSTEM  STATE  EQUATIONS 

The  system  equations  were  first  written  in  a  state  variable  format 
[6].   Thus,  in  matrix  notation,  the  system  state  equations  became 

X(t)  =  A  X(t)  +  B  U(t) 

Here  X(t)  was  the  matrix  of  system  state  variables  and  U(t)  was  the 

matrix  of  system  input.   From  the  overall  system  equations  (9-12),  Y(t), 

P„  (t) ,  X(t)  and  P  (t)  were  chosen  as  the  basic  problem  variables,  and 
Jta  m 

X   and  Fn    became  the  input  variables.   Thus,  the  applicable  state 
va      £  r 

variables  were 


Y  =  X(l)       X  =  X(4)       P  =  X(3)       P„   =  X(6) 

m  Jta 

Y  =  X(2)       X  =  X(5)       P  =  X(3)       P„   =  X(6) 

m  x,a 

Y  =  X(2)       X  =  X(5) 


U.  =  X         U   =  F0 
1    va       2    x. 


For  example  from  equation  (9)  one  can  write 

A 
X(2)  =FX(3)  -—  U  (13) 

In  a  similar  way,  the  system  equations  were  rewritten  in  matrix 
notation  (Figure  5) .   The  accompanying  signal  flow  graph  for  the  over- 
all system  equations  is  presented  in  Appendix  A.   The  state  equations 

were  solved  for  various  servovalve  stroke  values,  X   ,  and  for  various 

va 

loading,  F„,  using  the  IBM  simulation  language  CSMP  III  as  shown  in 
the  next  section. 
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V.   ANALYSIS  OF  SYSTEM  PERFORMANCE 

Manufacturer's  data  was  consulted  in  order  to  formulate  a  power 
piston  velocity  specification.   The  maximum  power  piston  velocity 
required  was  a  value  of  0.5  in/sec.   Common  operational  velocities 
ranged  between  0.4  -  0.5  in/sec. 

A.       PROCEDURE 

1.  Inputs 

Electro  hydraulic  servovalves  have  strokes  ranging  from  0.005 

-  0.010  in.  for  small  valves  up  to  10  gpm  capacity  to  0.015  to  0.030  in. 

for  50  gpm  valves.   In  this  case,  a  15  gpm  valve  was  used  and  input, 

X   ,  was  set  from  0.005  to  0.012  in. 
va 

An  external  loading,  FQ ,  was  assumed  to  range  from  0.0  to 
40,000.  lbf  in  the  simulations.   This  range  of  loading  could  be  expected 
in  operational  conditions. 

2.  Outputs 

The  scheduling  of  control  parameters  is  shown  in  Figure  6. 
Y  and  X  signified  the  displacement  of  the  power  piston  and  the  input 
piston  valve  rod  assembly.   P  was  the  load  pressure  at  the  power 
piston,  and  K_  was  the  fluid  spring  force  at  the  power  piston. 
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Figure  6.   Scheduling  of  Control  Parameters 


26 


B.   COMPUTER  SIMULATION  RESULT 

As  mentioned  previously,  computer  simulation  results  were  taken  for 

valve  stroke  displacements,  X   ,  from  0.005  to  0.012  in.   Figure  7 
r  va 

shows  the  resulting  power  piston  displacement.   As  can  be  seen,  the  sys- 
tem response  is  good. 

Figure  8  portrays  the  velocity  of  the  power  piston.   These  data 
demonstrate  that  the  velocity  of  the  power  piston  is  directly  propor- 
tional to  the  displacement  of  the  valve  input,  X   „   The  natural 

va 

hydraulic  frequency  of  the  system  was  approximately  95  Hz  for  the  con- 
trol unit  input  servovalve  and  piston  and  110  Hz  for  the  power  piston 
assembly.   These  calculations  are  shown  in  Appendix  B. 

Figure  9  is  a  plot  of  the  load  pressure,  P  ,  for  various  external 

loads.   Without  external  loading,  P  was  almost  zero.   After  an  added 

m 

external  load,  P  was  changed  slightly  for  a  range  of  Fp  values  from 
2000.0  to  40,000.0  lbf. 

Calculations  show  that  the  displacement  of  the  power  piston  upon 
loading  is  0.0016  in.  for  40,000.0  lbf.   Compared  to  the  maximum  move- 
ment (10.0  in.),  this  difference  is  about  0.016  percent.   Thus,  the 
displacement  of  the  power  piston  was  almost  independent  of  external 
loading. 

The  bulk  modulus  is  a  fluid  parameter  which  characterizes  the 
compressibility  of  the  liquid.   As  can  be  seen  from  Figure  10,  the 
value  of  the  bulk  modulus  has  a  slight  effect  on  power  piston  response 
time.   Between  bulk  modulus  values  of  100,000   and  300,000   PSI ,  the 
power  piston  response  was  changed  from  0.56  to  0.34  sec.   Compared  to 
a  maximum  required  time  of  19  to  20  sec,  during  which  the  power  piston 
moved  10.0  in.,  this  difference  in  power  piston  response  time  is  very 
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small.   This  is  an  encouraging  result  since  minimizing  entrapped  gas  to 
insure  high  3  values  is  a  difficult  and  costly  process. 

Figure  11  shows  the  influence  of  bulk  modulus  on  power  piston 
velocity  response.   Again,  in  view  of  a  total  system  response  time 
requirement  of  19  to  20  sec,  the  initial  oscillations  that  occur  within 
the  first  0.3  to  0.6  sec  of  system  operation  can  be  considered  negligible. 

Finally,  in  Appendix  C  is  shown  the  computer  program  developed  in 
this  study  and  a  typical  computer  output  (Figure  15) . 
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VI.   CONCLUSIONS 


In  this  study,  several  types  of  C.R.P.  systems  were  analyzed.   These 
included: 

1.  THE  TWO  VALVE-PISTON  COMBINATION  TYPE  WITH  SOLID  CONNECTING  ROD 
Results  for  this  system  have  been  presented  in  the  previous 
sections  of  this  report. 

2.  A  TWO  VALVE-PISTON  COMBINATION  TYPE  WITH  CONNECTING  FLUID  LINE  OF 

a.  circular  pipe  construction 

b.  annular  pipe  construction 

c.  double  annular  pipe  construction 

In  all  of  the  studies  involving  fluid  lines,  the  fluid  lines  had  very 
high  frequency  and  very  low  damping  ratio  characteristics.   Thus,  essen- 
tially the  influence  of  the  connecting  lines  on  response  could  be  neglec- 
ted.  Appendix  D  shows  the  result  of  various  pipe  line  response  calcula- 
tions . 

Computer  simulation  proved   to  be  a  useful  design  tool  regarding  the 
C.R.P.  system.   Full  scale  prototype  testing,  however,  is  the  only  con- 
clusive method  of  proving  C.R.P.  system  response.   Unfortunately,  testing 
is  both  costly  and  time  consuming.   Thus,  the  computer  simulation 
developed  in  this  study  could  be  used  to  designate  critical  testing 
instances  and  to  identify  those  areas  where  design  efforts  would  be  most 
productive.   In  this  way,  the  time  involved  in  optimizing  system  perform- 
ance as  well  as  subsequent  production  costs  could  be  reduced.   Likewise, 
the  present  computer  simulation  could  be  useful  in  specifying  the  correct 
approach  for  design  testing  practices. 
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APPENDIX  A:   SYSTEM  SIGNAL  FLOW  GRAPH 
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APPENDIX  B:   SYSTEM  RESPONSE 

The  governing  equations  could  be  separated  into  two  parts  -  the  input 

part  (servovalve  andpiston)  and  the  output  part  (power  piston) .   In  this 

Appendix,  the  hydraulic  damping  coefficient,  €,,    and  the  natural  frequency, 

CO,  /  of  each  part  are  determined, 
h 

A.   INPUT  PART  (SERVOVALVE  AND  PISTON) 

From  the  Figure  13,  the  characteristic  equation  is 

46  K         46  A  2 
C.E  =  s2  +  (_JL££)  s  +  (-JSLE-j 

V  M    V 

a  v   a 

2  2 

=  S   +  2£  CO,   S  +  CO, 

h      n 


where 


( 

2 

A 

46 

A 

fo, 

— 

V 

e   p 

h 

▼ 

M 

V 

V 
a 

K 

/6   M 

ca 

/    e   v 

t, 

A 

P 

1 

i^ 

Figure  13.   Signal  Flow  Graph  of  Input  Parts 
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For  nominal  values  of  bulk  modulus,  the  following  results 
pertain: 


3  (PSD 
e 

100,000 
300,000 


COh(Hz) 
95.77 
159.83 


0.00488 
0.00845 


B.   OUTPUT  PART  (POWER  PISTON) 

From  the  Figure  14,  the  characteristic  equation  is 


43  K      _     43  A  43  K  A 

c.e  =  s3  +  <— &-S)  s2  +  (-^-)  s  +  (— -S-SELJE.) 


v 


m 


mp   V 
x,   m 


v     Mo 
m     Z 


2  2 

=  (S  +  R)  (S   +  2C  w^S  +  a)   ) 

h     h 


Again,  for  nominal  values  of  bulk  modulus,  the  following  results 
were  determined: 


3  (PSD 
e 

100,000 
300,000 


0)  (Hz) 

110.45 
191.29 


0.00897 
0.04094 


Figure  14.   Signal  Flow  Graph  of  Output  Parts 


37 


APPENDIX  C:   COMPUTER  PROGRAM 


INITIAL 

CONSTANT  BE= 100000.0 

CONSTANT  FF=2000.0 

CONSTANT   CEM=1.5^ 

CONSTANT    QtfK=33^26.5^ 

CONSTANT   QAK=2057.08 

CONSTANT  XM=69.8?57 

CONSTANT  VM=4.9 

CONSTANT   V A  =7^9- 71 

CONSTANT  VC=l*4-326.0 

CONSTANT  AM= 1102.0 

CONSTANT  AA=57.67 

CONSTANT  W  =  l  1.938 

CONSTANT  PS =400.0 

CONSTANT  PSA =600.0 

CONSTANT  T1=0.0,T2=0.0,T3=0.0,T4=0.0,T5=0.0,T6=0.0 

DYNAMIC 

PARAMETER  Z=(0. 005,0. 01. 0.0 12,0.013) 

FC=FF*STEP(22.5; 

XIN1=Z*STEP(0.0, 

XIN2=2*STEP(5.0) 

XIN3=Z*STEP(10.0) 

XIN4=Z*3TEP(12.5) 
XIN5=Z*3TEP(15.0( 
XI N6=Z* STEP (20.0* 
XIN7=Z*STEP(25.0^ 
XIN8=Z*3TEP(30.0, 

x  va=xi  n 1-xi n2-xi n3+xi  n^-xi  n  5+xi  n6+xi  n7-xi n8 

ydd=(am/xm)*pm-( i./xm)*fc 

yd=intgrl(t1,ydd) 

y=intgrl(t2,yd) 

pmd=-(4.*be*qmk/vc)*y-(^.*be*am/vc)*yd-(4.*be*cem/vc)*pm  . . , 

+(^.*be*qmk/vc)*x 

pm=intgrl(t3,pmd) 

fk=0.86*w*(ps-pm)*(x-y) 

xdd=(fk/vm)*y-(fk/vm)*x+(aa/vm)*pla 

xd=intgrl(t4,xdd) 

x=intgrl(t5,xd) 

plad=-(4.*be*aa/va)*xd-(^.*be*cea/va)*pla+(^.*be*qak/va)*xva 

pla=intgrl(t65plad) 

OUTPUT  Y 

PAGE   MERGE 

METHOD  RKSFX 

TI MER  FI NTI M=35 . 0 , OUTDEL=0 . 5 , PRDEL=0 . 5 , DELT=0 . 002 

END 

STOP 
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APPENDIX  D:   FLUID  LINE  RESPONSE 

When  a  fluid  is  used  for  transmitting  pressure  variations,  it  is 
necessary  to  determine  the  response  of  the  transmission  system.   In 
the  course  of  this  study,  the  following  line  response  determinations 
were  formulated. 
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Figure  16.   Circular  Type  of  Line  Response  Determinations. 
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Figure  17.   Annular  Type  of  Line  Response  Determinations. 
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